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ABSTRACT 



The effects of gas pressure gradients on the motion of sohd grains in the solar nebula substantially 
enhances the efficiency of forming protoplanetary cores in the standard core accretion model in 'hybrid' 
scenarios for gas/ice giant planet formation. Such a scenario is enhanced core accretion which results 
from Epstein-drag induced inward radial migration of mm-sized grains and subsequent particle subdisk 
gravitational instability needed to build up a population of 1 km planetesimals. Solid/gas ratios can 
be enhanced by nearly lOx over those in Minimum Mass Solar Nebula (MMSN) in the outer solar 
nebula (a > 20 AU), increasing the oligarchic core masses and decreasing formation timescales for 
protoplanetary cores. A 10 core can form on ^ 10^ — 10^ year timescales at 15 - 25 AU compared 
to ^ 10® years in the standard model, alleviating the major problem plaguing the core accretion model 
for gas/ice giant planet formation. 

Subject headings: solar system; formation planetary systems: formation planetary systems: proto- 
planetary disks 



1. INTRODUCTION 

The main challenges facing the standard core accre- 
tion model for planet formation are finding a mechanism 
for building up a population of planetesimals and form- 
ing sufficiently massive proto-Uranus and proto- Neptune 
cores before the end of the oligarchic growth stage and 
before the solar nebula dissipates on 10^ yr timescales. 

In the core accretion model, planetesimals are built up 
from approximately /im-sized grains to km-sized bodies 
by coUisional sticking. Two major problems face col- 
lisionally sticking these bodies together. First, grow- 
ing particles from centimeter through kilometer sizes is 
problematic. While many have argued for the coagula- 
tion of grains based on some requisite 'coagulation veloc- 
ity' (e.g. Weidenschilling 1997) such velocities are often 
much higher than those below which sticking might oc- 
cur based on microgravity experiments such as Blum & 
Muench (1993). More recently, Colwell (2003) showed 
that at least in some conditions cm-sized grain coagula- 
tion can occur, but this occurs for relative velocities be- 
low Ylcmj s and greater relative velocities are encoun- 
tered in either a turbulent or laminar disk (see Youdin 
& Shu 2002, hereafter YS02, and references therein). An 
operative sticking mechanism also remains in question: 
both known solid state sticking forces and the bodies' 
self gravity are arguably too weak (YS02). Second, even 
if such sticking were possible, the timescale for meter- 
sized bodies to spiral in to the Sun (by gas drag) is 
short (Weidenschilling 1977), 10^ years at 1 AU, such 
that planetesimal formation in all but the outermost re- 
gions of the solar nebula might not proceed fast enough. 
One can circumvent this problem if a particle subdisk 
gravitational instability (GI) can very quickly build up 
a population of Ifcm planetesimals directly (Goldreich & 
Ward, 1973; YS02). 

However, formation timescale problems still remain for 
ice giants Uranus and Neptune. For a MMSN distribu- 
tion, the oligarchic growth stage does not yield Neptune- 
Electronic address: lcurrie@astro.ucla.edul 



mass cores (for a > 15 AU), so that a 15M^ core must 
be formed by colliding together sub Earth-mass oligarchs. 
Levison & Stewart (2001) showed that the embryos' mu- 
tual perturbations result in a large number of ejections 
(not mergers) of embryos in the outer solar system, so 
that the formation timescale for a 10 — 15M0 proto- 
Neptune core is prohibitively long. The two most widely 
explored attempts at a solution are positing a massive 
(> 6x MMSN) disk during the planet formation epoch 
(Goldreich et al., 2004, hereafter GO04) or that ice giant 
cores were gravitationally scattered to their present or- 
bits after brisker formation in the trans-Saturnian region 
(Thommes et al. 2002; Thommes et al. 2003). However, 
these ideas, while plausible, introduce very case-specific 
assumptions about the initial state of the solar nebula or 
its dynamical history. 

The purpose of this paper is to show that another 
mechanism for forming ice giant cores within the core 
accretion model exists that is more generally applicable 
to protoplanetary disks. Specifically, one may be able to 
form Uranus and Neptune-mass cores in situ and on solar 
nebula dissipation timescales from what is currently the 
most plausible mechanism for inducing particle subdisk 
GI necessary for planetesimal formation. In Section 2 
& 3 we first review this planetesimal formation mecha- 
nism, called here the Sekiya- Youdin gravitational insta- 
bility model, and then argue that it could aid ice giant 
formation in proceeding quickly. Next, in Section 4 we 
describe our recalculation of grain pileups in this model 
originally done by YS02 and Youdin & Chiang (hereafter 
YC04) , develop a simple model for calculating protoplan- 
etary core masses at the end of oligarchy and core forma- 
tion timescales. In Section 5 we present our the results 
of our core mass and formation timescale calculations, 
showing how the Ep, solid body surface density, profile 
resulting from this planetesimal formation mechanism al- 
ters these values substantially. In the discussion section 
we put this model in a historical context, summarize it, 
and describe what research within this model should be 
done further. Finally, we include an appendix describing 
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the protoplanetary disk conditions necessary to allow GI. 

2. PLANETESIMAL FORMATION FROM SEKIYA-YOUDIN 

GI 

Unless one assumes a hitherto unknown efScient stick- 
ing mechanism for growing 1cm — 1km bodies one likely 
has to posit some sort of particle layer GI to build up 
km-sizcd planetcsimals. Planetcsimal formation by GI, 
as developed by Goldreich & Ward (1973), was long con- 
sidered problematic, though. The chief reason for this 
was that as particles settle towards the disk midplane the 
disk becomes sufficiently vertically stratified that Kclvin- 
Helmholtz turbulence develops. Turbulent eddies result- 
ing from the Kelvin-Helmholtz instability induce parti- 
cle random velocities that are too high to allow GI, and 
equilibrium vertical profiles for particles result in spatial 
densities that are at least an order of magnitude too low 
(Weidenschilling 1980; Cuzzi, Dobrovolskis, & Champ- 
ney 1993 ; Weidenschilhng 1995). 

Sekiya et al. (1998), however, found that if the ratio 
of dust to gas surface densities was sufficiently enhanced 
over standard solar values, turbulence induced by Kelvin- 
Helmholtz vertical shear will not be sufficient to stir all 
the solids: the rest will be able to precipitate towards the 
disk midplane and induce GI. YS02 and YC04 provide 
a likely mechanism for achieving this enhancement by 
pileup of grains from gas drag-induced migration, avoid- 
ing the main difficulty with GI from Goldreich & Ward 
(1973). We now describe this mechanism following YS02 
and YC04. 

Because gas is more sensitive to pressure gradients, 
it orbits at slightly sub-Keplerian velocities where the 
relative velocity difference between gas and dust is pro- 
portional to 1] = —{dPg / dlna) /2pgvf. ^ (cg/vk)"^. The 
dust then experiences a headwind, and thus a drag 
force. For a stopping time of a particle moving rel- 
ative to the gas of Tgtop = Pp^/pg'^g^ ^^'^ S^^ ^ind 
particle densities of pg and pp respectively and parti- 
cle sizes s, the inward particle flux through the disk 
is / = fep + fturbi migration due to Epstein drag is 
/ep = Pp{^fvep,ind = Pp{Pg / p)'^'^V^TstopVk and turbu- 
lent stresses is fturb ~ PpVep,ind ^-§^iPpz/p) (whose ef- 
fect is small: see YC04). This inward particle flux from 
Epstein drag-induced migration results in grain pileups if 
the flux decreases with decreasing stclloccntric distance 
a. Explicitly, whether such a condition is satisfied de- 
pends on the radial profiles for solid and gas surface den- 
sities Sp and T,g as well as the temperature profile. 

In a MMSN model, one starts with solid body gas pro- 
files of Sp oc a"" and Eg oc a^P and a temperature profile 
of T oc a~^. The migration rate due to Epstein drag is 

Va ~ Vep,ind = 2^r]Tstop^'^a. (1) 

As Va depends on rj and tgtop which in turn depend on 
T and S^. Isolating the radial power law dependence of 
T and Eg on a one finds that the drift rate dependence 
goes as Va oc a'^, where 

d = p-q + l/2. (2) 

The mass accretion rate is oc EpOWa oc where E = d — 
n+1 01 E — ^+p—n—q for grain migration in the Epstein 
drag regime . If E > initially then the accretion rate 



increases with stclloccntric distance, resulting in particle 
pileups and thus Ep enhancement. Epstein drag then 
results in grain pileups from the particle flux's sensitivity 
to the gas density and temperature profiles. 

Once enhancements of Ep/E^ from migration are suf- 
ficiently high such that the solid particle surface density, 
Ep, rises above some critical value Epc then particle sub- 
disk GI (and thus planetcsimal formation) commences 
from dust in excess of Epc. Epstein drag- induced migra- 
tion brings about particle layer GI because pileup of solid 
grains from migration can raise Ep above Epc, overcom- 
ing Kelvin-Helmholtz shear instabilities. Particle layer 
gravitational instabilities are allowed once the solid to 
gas ratio is enhanced by 5 — 20x (see YC04). Garaud 
& Lin (2004) find a similar criteria for GI using a more 
sophisticated two-fluid treatment. 

3. ENHANCED EFFICIENCY OF PROTOPLANET CORE 

ACCRETION FROM SEKIYA-YOUDIN GI 

YS02 and YC04 showed how Epstein drag-induced 
grain migration can procure planetesimal formation but 
didn't explore the effect that such migration has on the 
efficiency of planet formation. Epstein drag-induced mi- 
gration of grains, as shown by YS02 & YC04, yields 
substantially different solid body surface density profiles 
from an MMSN profile, enhancing Ep. Since the cri- 
teria for GI (and planetesimal formation) is a 5 — 20a; 
enhancement of Ep, the mass of planetcsimals available 
form which to form protoplancts should greatly increase 
within the dust subdisk as the mass from which such 
planetesimals are formed is locally enhanced by import- 
ing material from the outer solar nebula via migration. 

As the masses of post-oligarchic cores as well as pro- 
toplanet formation and inelastic coUisional damping and 
collision timescales depend strongly on Ep a seemingly 
simple enhancement of Ep by about an order of mag- 
nitude could have very important implications for the 
efficiency of planet formation in regions with such en- 
hancement. The effect wasn't obvious from YS02 and 
YC04 since both papers did not include solid enhance- 
ment from water ice condensation in their initial Ep pro- 
files. For example, model H in YC04 then triggers GI 
at an outer boundary of only ~ 2 AU while model Af 
triggers GI at an inner boundary with an enhancement 
of only ~ 1.5 — 2 X T,p^MMSN at 15-25 AU when one 
compares the resulting Ep profile with an MMSN pro- 
file including the water ice enhancement. However, if 
one includes water ice condensation in Y^p^init then pileup 
should proceed faster further away from the Sun. Plan- 
etesimal formation could then be triggered earlier and 
the resulting Ep profile could influence planet formation 
efficiency in the outer solar system (10-30 AU). 

4. CALCULATIONS 

We then are motivated to redo the YS02 calculations 
(equations 22-30) for the radial migration of mm-sized 
grains but include the solid enhancement for a > 3AU 
due to the water ice condensation. We also estimate core 

masses expected at the end of the oligarchic growth stage 
as well as formation timescales for lOM^ cores. 

4.1. Recalculation of Grain Radial Drift/Pileup 

Following YS02 we assume a population of uniformly- 
sized particles comprising the solid body column Ep. Its 
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time-dependent profile is 



9Ep dY.p T.p d 



(3) 



which has an analytical solution: 

Sp(a,i) = 5]„a-i-'^af+i-"(a,t) (4) 

where d = l/2p — q and ai{a,t) is the initial location of 
a particle now at a after time given by Eq. 26 of YS02: 

a 



ai{a,t) = a(l - {d - 



(5) 



4.2. Oligarchic Core Masses and Core Formation 

Timescales 

We then estimate oligarchic core masses and core for- 
mation timescales in the standard MMSN model and 
the model modified for migration of mm-sizcd grains. 
Oligarchic growth ends in once oligarchs' velocity ex- 
citations from mutual viscous stirring overtakes their 
damping rates which result from dynamical friction with 
smaller planctcsimals. This condition is met once the 
"surface density" of protoplanets, Sm = Mpro/'^.TraSa 
where 5a is the feeding zone size, becomes comparable to 
the residual surface density of smaller bodies such that 
Sjvf ~ ^m{t) (GO04). Assuming that little mass is yet 
lost from the system this condition can be related to the 
initial surface density of solids: 

EM + S„(t) = i;p (6) 

Solving for Mp^o then yields the protoplanet mass at the 
end of ohgarchy {MoUg): 

Moiig ~ Qm{b/l{)f-^{Y.p/lQgcm-^f-^{a/lAUfM^. 

(7) 

After oligarchy ends dynamical interactions between oli- 
garchs results in ejection of many of them, not mergers 
(Levison & Stewart 2001). 

For calculating the core formation timescale we gener- 
ally follow Ida & Lin (2004) and Kokubo & Ida (2002) 
with some modifications depending on the size of ac- 
creted planetesimals gleaned from Rafikov (2004). From 
Kokubo & Ida (2002), the mass accretion rate of a pro- 
toplanetary core depends on Sp, the mass of the core at 
time t, the core's size Vp, its stellocentric distance a, and 
the velocity dispersion of accreted planetesimals: 



M = CttS, 



2GMr„ 



P 9 i 9/-.' 



(8) 



where < >5 is the rms eccentricity of planetesimals 
and C is a factor of order unity. The planetesimals' rms 
eccentricity is found by c^iTiating the; viscous stirring and 
gas drag timescales Tyg — Tgas (Kokubo & Ida (2002)). 
The formation timescale is then given as Tc^acc = ^ • The 
approximate formation timescale for a core is then given 
by Ida & Lin (2004): 

Tc,acc 1.2 X 10^(105cm"VSp)(2400,gcm-2/Sg)"-'' 

x(a/lAC/)°-<^(Me/Me)i/^(m/10i*'ff)2/i52/ears,(9) 

where m is the mean accreted planetesimal mass for 
a protoplanet in a core feeding zone size of IQRhui 
{b = 10). The rate of accretion of ~ lOOm — lOfcm 
planetesimals can be described by the above equation 



and accretion is then said to be 'dispersion dominated'. 
However, if the accreted planetesimals are small enough 
their relative velocities to cores are set by differential 
shear within the Hill radius of the core. In this case, the 
approach velocity of a planetesimal to the core is ~ QRh 
where Rh is the core's Hill radius, and the vertical com- 
ponent of the planetesimal's velocity, v^, small. For a 
small enough Vz the planetesimal disk becomes very thin 
and the embryo can accrete the entire vertical column 
of planetesimals (Rafikov 2004). This regime of accre- 
tion is called 'shear- dominated accretion' as the shear 
from the planetesimal disk around a core sets that core's 
velocity. Rafikov (2004) then shows that the formation 
timescale for protoplanetary cores accreting shear- dom- 
inated planetesimals then is shorter. We rewrite the 
timescale here, generalizing it to variables Sp, pp, M^, 
and X- 

T,,acc 2.36 X 10^x'\lOgcm-yj:p) 
x{M,/M^Y/\pp/lgcm-^Y/^years, (10) 

where % is the fraction of planetesimals which are shear 

dominated and pp is their mass density. These isolated 
body equations for MoUg and Tc,acc should be valid as 
long as growth occurs prior to the end of oligarchy and 
on timescales comparable to or less than the complete 
dissipation of the gas disk. 

4.3. Size of Accreted Planetesimals 

The size of accreted planetesimals is likely to be smaller 
than the ~ 1km bodies formed after GI. Kilometer- 
sized objects inelastically collide and fragment if their 
collisional velocities exceed Qd, the energy per gram 
needed to release half of the planetesimal's mass upon 
collision. For a 300 m projectile hitting a 1km target 
~ 4 X IQ'^erg/g (Leinhardt & Richardson 2002), 
where Qd decreases for mass ratios approaching unity: 
low Qd's are consistent with the standard low internal 
strength, 'rubble pile' model for icy fcm-sizcd bodies such 
as comets (e.g. Asphaug & Benz 1996). The random ve- 
locity needed to disrupt a 1km planetesimal is ~ lOm/s 
for a comparably-sized projectile of s = 500m. This 
velocity is comparable to vh at 30 AU. The velocity 
of collisions actually experienced between ^ l/cm-sized 
planetesimals near a protoplanetary core is set by the 
Hill velocity, vh, because their random velocities are set 
by scattering encounters with the core (GO04, Kenyon 
& Bromley 2004, and Rafikov 2004). Furthermore, ran- 
dom velocities after a scattering encounter are likely to 
exceed vh (see Rafikov 2004, equations 23, 27, 42, and 
46) unless the bodies already have s « 1km. This 
means that 1 km planetesimals scattered by protoplan- 
etary cores should fragment upon collision with one an- 
other, an outcome confirmed by numerical simulations of 
planetesimal interactions around cores (e.g. Kenyon & 
Bromley 2004). 

This is important because, as was alluded to before, 
the size of planetesimals very fundamentally affects core 
accretion rates. Specifically, for small planetesimals ac- 
cretion can be much more rapid when planetesimal ac- 
cretion occurs in the aforementioned 'shear-dominated' 
regime. We now describe the conditions for such accre- 
tion as they are given in Rafikov (2004). 

Whether planetesimal accretion is shear or disper- 
sion dominated depends largely on values for a fiducial 
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mass, Mf and a fiducial size, rf which relate the plan- 
etesimal mass and sizes to planctesimal velocities in- 
duced after embryo-planetesimal scattering events. In 
particular, Mf = Mqtj^ ~ 8 x 10'^A{a/lAUf/'^g and 



0.2 



Ps C5 



120(a/ly4;7) 



-7/4. 



Also important 



is a size, Tg, where the Reynolds number of the gas, 
Re, is Re = Ret = 20, corresponding to changes in 
the drag coefficient experienced by a planctesimal mov- 
ing through gas and where the planetesimal's velocity 
dispersion equals its velocity difference with the gas: 
V = 6vg(Ra,&kov 2004). This size is given as 

n = ^^^2.5{a/lAUf/'m, (11) 

3 Cs 

where A ^ 1.25{a/lAUy^^^cm: the molecular mean free 
path (Sg is reduced by a factor of 1.25 compared to 
Rafikov (2004) affecting r/, A, and r^). 

As shown in equations (42) and (46) of Rafikov (2004), 
the velocity of a planctesimal relative to the Hill velocity, 
vh, depends on the core mass, Afc, as well as on Mf, A, 
rf, and rg. li v < Vh then shear-dominated accretion 
sets in, and thus we can describe the boundary between 
dispersion and shear-dominated accretion as it depends 
on these parameters. For the Stokes drag regime the 
condition for shear-dominated accretion is 



Tf S 

and for the Epstein drag regime it is 

A s 



(12) 



(13) 



Thus, after some time of planctesimal accretion, the 
planetesimals must collide frequently enough such that 
the population is ground down, keeping s small, so that 
shear dominated accretion may persist. 

4.4. Model Assvmptions and Initial Conditions 

For the radial migration calculations we begin with 
the following distributions of solids and gas. We set 
So = e X 10gcm~^ where e = 1 interior to 3AU, and 
e = 4.2 exterior. The initial gas & solid surface den- 
sities arc given as Eg = 2400 x (a/lAU)~Pgcm~'^ and 
Ep = So X (a/lAU)^" and temperature profile given as 
T = 280 X {a/lAU)-m where p & n = 1.5, and q = 0.5, 
corresponding to Model H in YS02. The disk initially 
extends out to 250 AU in this model, though the outer 
edge of the solid component of the disk moves inward 
with time. For this model, the drift rate can be given as 



3{a/lAU) 



15 IQppS lAU 
gcm~'^ lO^yr' 



(14) 



We assume that the grains are comparable in size to the 
largest chondrules (.s = 1mm) to make direct compar- 
isons for the same nm in YS02. We do calculations for 
the formation timoscalo of a lOMg core, one with the 
mean size of accreted planetesimals, s, set at 1 km, an- 
other with s = 100m to take into account fragmentation, 
and s = 10m to further explore shear-dominated accre- 
tion. The density, pp, is set to 'Sgcm^^, representing an 
upper limit for planetesimals and ice-rock cores. 



5. RESULTS 

The solid particle migration run was stopped once one 

part of the disk had the requisite enhancement to in- 
duce a particle layer gravitational instability according 
to YS02. Figure 1 shows the time evolution of Sp(a,i) 
from an initial MMSN profile to the onset of particle sub- 
disk GI. When the bump from water ice condensation is 
included in the surface density profile the time for the 
onset of GI is ~ 10 x shorter and occurs much further 
out at ^ 27 AU than in the same model in YS02 and 
YC04, yielding Sp(a, t)/Sp(a, 0) ~ 6 - 8x from 20-25 
AU and ~ 9.15x at the outer edge of the particle sub- 
disk {26.7AU) before GI is initiated at t = 4.25 x 10^ yr. 
In the 'enhanced core accretion' scenario, the disk con- 
tains far more material from which to form protoplanets 
in precisely the regions where such formation has been 
most problematic for the standard core accretion model. 
As the figure shows, this occurs because the initial distri- 
bution in Sp is now 'squeezed' into a smaller surface area 
bounded roughly by the ice giant planet forming region 
due to Epstein drag-induced migration. 

Figure 2 shows the masses of oligarchs for standard 
and migration-enhanced core accretion models (labeled 
Sekiya-Youdin migration). For the latter, the oligarchic 
mass for Jupiter at 5.2 AU is 3.4^^®, 21.1Af0 for Uranus 
at 20 AU and BTMg, for Neptune at 25 AU, values 
much greater than in the standard model, implying that 
Uranus and Neptune could have accreted most, if not 
all, of their mass (15 — 17M0)by the end of the oli- 
garchic growth phase. The masses of embryos at the 
end of oligarchy are large enough that additional, post- 
oligarchic growth can be negligible. Dynamical fric- 
tion exerted by smaller bodies during oligarchic growth 
should dampen the embryos orbital eccentricities inhibit- 
ing ejection. Substantial post-oligarchic growth by co- 
alescing is then unnecessary since cores have already 
grown sufficiently massive during oligarchy. Extra proto- 
planetary mass in the outer solar nebula is ejected after 
oligarchy ends (GO04). 

Figure 3 compares formation timescales for a lOM® 
protoplanct core for the standard core accretion model 
and the model utilizing the particle subdisk GI from 
Sekiya and Youdin's work to set Sp. For the standard 
MMSN surface density profile the formation timescale for 
a lOM© proto-Neptune core is 4.4 x 10^ and 1.75 x 10* 
years at 25 AU for s — 1km and s — 100m respec- 
tively. Since Tc,acc>Tdissip and MoUg is significantly less 
than Mneptune, the MMSN model will likely fail to pro- 
duce Neptune-like planets on requisite timescales. How- 
ever, for the enhanced model the formation timescales 
arc 5.5 x 10'' and 2.2 x 10*" years at 25 AU for s = 1 km 
and s = 100m. It then could be possible to form proto- 
Neptune nearly in situ on a timescale less than the disk 
dissipation timescale (10-30 Myr) for s < 100m. 

Furthermore, a formation timescale of 20 Myr at 25 
AU is most likely an overestimate as, using arguments 
from Rafikov (2004), protoplanets here should grow a 
large fraction of their mass while the accreted planetes- 
imals are in the fast, shear- dominated regime. Initially 
one starts with 1-10 km planetesimals where the initial 
growth rate is fast. However, once the core mass in- 
creases enough, planctesimal accretion becomes disper- 
sion dominated and growth occurs much more slowly. 
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To get back to faster, shear-dominated accretion, the ini- 
tial population of planctesimals must drop in size enough 
such that s « 1km. Quantitatively, we refer back to 
the equations (11) through (13). At 25 AU M/ and r/, 
are ~ lO^^g and 0.43to respectively and rg is ~ l.Skm. 
Then, in the Stokes regime, dispersion-dominated accre- 
tion persists for Me > 3.8 x lO^^g and 3.8 x lO^^^f for 
s = 1km and 100m the latter still not an appreciable 
fraction of Neptune's current mass. 

However, the ^ 7 — 9.x enhancement resulting from mi- 
gration and pileup of mm-sizcd grains and subsequent GI 
helps to bring one back into the shear-dominated accre- 
tion regime earlier and persist there for longer times by 
the following argument. This enhancement drives down 
the collision timescale by nearly an order of magnitude as 
the collision rate is proportional to the number density 
of planctesimals and thus Sp (e.g. Wetherhill & Stew- 
art 1993, equation Al). This increases the dominance of 
collisions and collisional damping over gas damping and 
inspiral of solids (the latter timescale is longer than gas 

damping by a factor of 1/(1. 6x lO^'^aaM^) ~ 125). IfSpis 
enhanced we should expect more frequent collisions early 
on. As the collisions tend to be disruptive at expected 
impact speeds of f > w_f/, an enhancement in Ep should 
result in planetesimal swarms of smaller mean sizes. 

The smaller sizes translate into planetesimal dynam- 
ics occurring in different drag regimes and longer times 
spent in shear dominated planetesimal accretion. Specif- 
ically, at 25 AU, planctesimals are in the Epstein regime 
for s < A ^ 88m. As this is slightly below an order 
of magnitude reduction in the size of the initial swarm 
and collision rates between planctesimals are greater 
for enhanced Ep, it is likely that an appreciable frac- 
tion of accreted solids are in the Epstein drag regime 
and this regime is reached earlier than when a stan- 
dard MMSN Sp distribution is used. One stays in 
the shear-dominated, Epstein drag regime as long as 
Me < 5 X lif^'g and 6 x 10"^^ g for s = 50m and 10m, 
the latter being c;omparable to the present rock/ice mass 
of Neptune. In the context of planetesimal accretion this 
means that cores can accrete planctesimals in the fast 
shear-dominated regime through a larger fraction of their 
total mass. This has a drastic effect on the core forma- 
tion timescale. For shear-dominated accretion, we then 
find the formation timescale for a lOM® core at 25 AU 
to be T ~ 2.26 x lO'^x^^yr = 2.26 x lO^yr if 10% of plan- 
ctesimals are shear dominated (shown in Figure 3) and 
2.26 X lO^yr if all of them are: both values are less than 
the solar nebula dissipation timescale, T^issip ~ lO^yr. 
Formation timescalcs in the outer solar nebula are then 
sufficiently reduced and conditions are more likely to 
be met such that accelerated protoplanet growth per- 
sists through a sufficient fraction of proto-Neptune's final 
mass. 

6. DISCUSSION 

The general idea in this paper that an enhancement of 
solid material in the Uranus & Neptmie regions of the 
solar nebula during the first 10 Myr is what led to rapid 
formation of ice giants is not new. Various incarnations 
of this idea have appeared in the last twenty years (e.g. 
Lissauer 1987; Pollack et al.l996; GO04) and even when 
such an enhancement was not considered the primary 
factor in forming Neptune on gas dissipation timescales. 



Tdissip, it often plays a supporting role in increasing a 
proto-Neptime's feeding zone zone (e.g. Bryden et al. 
2000). These papers, however, usually assume such an 
enhancement of solids over MMSN values results solely 
from the disk mass itself being much greater than MMSN 
values. This paper relaxes that requirement, assuming a 
Eg and initial Ep profile comparable to MMSN values 
(total disk mass here is ^ O.O5M0 spread over 250 AU, 
out to 30 AU it is ^ O.O2M0), with very large enhance- 
ments of the solids achieved from grain migration and 
GI. Thus, in this model it may be possible to get the 
same benefit you have from an initially massive disk (e.g. 
more solids from which to form planets, higher planetes- 
imal collision rates, etc.) without the 'cost' of having 
that disk mass be unreasonably large. 

That a global radial redistribution of solids influences 
some part of the planet formation process has also been 
suggested previously. Specifically, Stepinski & Valageas 
(1997) and Kornet et al. (2001) posited that gas drag in 
a turbulent disk result in solid surface density distribu- 
tions that differ from the gas surface density, providing 
enhancements in the inner regions of the solar nebula 
when compared to their initial distributions. These pa- 
pers, though, use a very simple treatment for grain coag- 
ulation for 1cm — lOfcm bodies, the size regime where a 
hierarchical sticking hypothesis seems to be highly prob- 
lematic. The effect any solid body redistribution has 
on the efficiency of post-planetesimal growth wasn't ex- 
plored: rather, at least for Kornet et al. (2001) emphasis 
was placed on how different initial conditions for the disk 
mass and radial extent explain the diversity of planetary 
systems via their differences in grain redistribution. It 
appears that grain redistribution can lead to not only a 
diversity of system architectures but changes in the effi- 
ciency of planet formation. 

The way in which planetesimal formation, particu- 
larly that of the GI model for planetesimal forma- 
tion. influences the efficiency of planet formation has been 
hitherto underexplored and warrants further investiga- 
tion. A comparison between analytical estimates of core 
formation timescales and numerical calculations assum- 
ing a Ep distribution set by dust migration models such 
as YS02 should be made. Furthermore, incorporating 
the 'shear-dominated' planetesimal accretion model into 
such calculations is needed as many papers making pro- 
nouncements about the efficiency of protoplanetary core 
formation (e.g. Thommes et al. 2003) adopt a pre- 
scription for the mass accretion rate similar to that for 
'dispersion-dominated accretion' in this paper, where a 
planetesimal's rms equilibrium eccentricity is set by bal- 
ancing viscous stirring by the protoplanet with gas drag 
from the solar nebula. As argued here and especially in 
Rafikov (2004), if the planctesimals are small then accre- 
tion might proceed more rapidly. 

If one considers the effect that Epstein drag-induced 
migration has on the distribution of solids, the forma- 
tion of cores in the outer solar nebula occurs on far 
shorter < lO^'yr timescales than in the standard core 
accretion model. Assuming an initially massive disk or 
requiring substantial gravitational scattering of proto- 
Uranus and proto-Neptune cores, while still plausible, 
may be unnecessary for forming Uranus and Neptune in 
< lO^years. Then 'hybrid' mechanisms such as this one 
resulting from Sekiya-Youdin particle subdisk GI could 
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eliminate core formation timescale problems in the outer 
solar system. Furthermore, the mechanism outlined in 
this scenario is generic if planetesimal formation proceeds 
by particle subdisk gravitational instability as described 
by YS02 & YC04, operating regardless of whether oth- 
ers do. Gas pressure gradients then play an indispensible 
role in forming the cores of gas/ice giant planets quickly, 
as they do in the rival disk instability model (see Boss 
et al. 2002, Haghighipour & Boss 2003), though in a 



slightly more indirect (and less obvious) way. 



The author would like to thank Andrew Youdin, Scott 
Kenyon, and Richard Nelson for extremely fruitful dis- 
cussions and Nader Haghighipour, Brad Hansen, and the 
anonymous referee for suggestions that strengthened the 
arguments presented in this paper. This research was 
supported in part by the NASA Astrobiology Institute. 



APPENDIX 

DISK CONDITIONS FOR GRAIN SETTLING/GI 

It is important to ask whether the conditions for particle subdisk gravitational instability can be met. Specifically, 
the main requirement for GI is that the disk must be sufficiently passive to initially allow the dust layer to collapse 
to a thickness set by a balance between settling and Kelvin- Helmholtz shear. If viscous stirring is strong enough to 
prevent this amount of settling then the entire GI mechanism for forming planetesimals becomes highly problematic. 
We show here that the assumption that GI can occur is at least reasonable. 

We first analyze previous attempts at constraining the equilibrium dust layer scale height achieved by balancing set- 
tling with diffusivity generated by Richardson turbulence. YC04 calculated this criteria to be a < 10^^(a/lj4t/)'^+P~* 
or a < 10^''{a/lAU)^ for a MMSN model. The coefficient a must then be less than ~ 6.3 x 10^'^ for a — 25 AU. More 
stringent requirements can be found from Dubrulle et al. (1995). From their equations 36 and 37, an equilibrium 
distribution of solids is written in terms of the ratio between the dust and gas scale heights (h and H), the friction 
(stopping) time Tstop, and a: 

~ ^Ja/nTstop■ (1) 

The criteria for GI in terms of volume density is pcr ~ Mq/q^ or ~ 2.4:XlO~^^ gcm~^ or about a factor of 1000 x increase 
over its solar abundance (and ^ lOOx the gas density). We then can rewrite the above condition as (S^/lOOSg)^ ~ 
10~^ ~ a/{Tstop^)- As Tstop^ ~ 10~*(a/lA{/)^'^ the criteria then becomes a < 10~^ for grain sedimentation to a 
requisitely thin dust subdisk. While estimates for a from molecular viscosity are small enough 10^^^), values for 
a needed to explain mass accretion rates onto T Tauri stars range from 10^^ — 10~^. 

Thus at first glance it seems as though conditions for GI include an implausibly low disk viscosity. However, as 
pointed out by YC04, turbulence need not be isotropic, and when a is resolved into components the criteria radially is 
given as a,. < ^fsiop ^ 10~^{a/lAU)'^^'^ (from their equation B2). The condition at 25 AU is then that < 1.25 x 10~^ 
which is reasonable. The effective a,, may differ substantially from the if the viscosity source is more effective at 
mixing radially than vertically. An example of a ratio between horizontal and vertical viscosities of order 10*^ is cited 
by YC04 for Earth's oceans and atmosphere. 

A stronger argument can probably be made based on investigating disk structure models in terms of whether a 
particular region is or is not viscously evolving, especially considering the effect of dust grains on levels of disk viscosity. 
The disk may spend a substantial fraction of its lifetime accreting through only thin surface layers sandwiched between 
a quiescent, magnetically dead layer containing a large fraction of the total column of the disk material in so called 
'layered disk' models (Gammie 1996). The anomalous viscosity then presumably tapers off and disappears beyond a 
certain gas column density into the disk a result consistent with MHD simulations of layered disks from Fleming & 
Stone (2003). Fromang et al. (2002) shows that for a = 10~^ and M ^ 10~^ Mq / yr the outer radius of the dead zone 
can be as large as 100 AU and occupy 90 % of the vertical column density of the disk for metal fractions, xm, slightly 
below cosmic abundances. For a = 10~^ a dead zone exists for all values of the disk metal fraction, not just subsolar 
values. Furthermore, these calculations were done assuming that the grains do not affect the magnetic coupling to 
the disk, and the authors specifically assumed that settling of the grains to the disk midplane had already occurred: 
thus their ability to inhibit coupling of magnetic field lines to the gas would be irrelevant except for a vanishingly thin 
dust layer at the midplane. Fleming & Stone (2003) also suggest that including the effects of grains would inhibit the 
ability of magnetic field lines to couple to the gas. 

In absence of this assumption, Fromang et al. (2002) suggests that because grains are effective at 'scavenging' charge 
they may induce a magnetically dead state within the subdisk they are contained, preventing turbulence from the 
anomalous source of viscosity from stirring them up and thus allowing sedimentation to the disk midplane if vertical 
mixing isn't particularly strong. This intuition now seems to have been confirmed numerically. Recently, Nelson et al. 
(2005, in prep) has investigated the effect that dust grains have on the coupling of magnetic field lines to the gas in the 
disk: thus their effect on the level of MRI turbulence. This work suggests that dust grains are particularly efficient at 
scavenging charge in the disk, and their inclusion leads to dead zone Sg's that are much larger than those in Fromang 
et al. (2002). They find that in order to drive MRI turbulence throughout the disk cross section at some distance, one 
needs to reduce the number density of dust grains by a factor of 10*. Accretion through the entire vertical column of 
the disk would then seem to presume that substantial grain growth and settling had already occurred. Clearly, one 
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Fig. 1. — Evolution of the solid body surface density profile, Ep. 
For this model (equivalent to Model H in YS02 & YC04), Epstein 
drag experienced by mm-sized grains causes them to migrate, pile 
up, and induce particle subdisk GI to form planetesimals. This 
leaves a ~ 8-9x enhancement in the outer regions of the particle 
disk after ~ 4.25 X 10^ yr, yielding a Sp diflfering significantly from 
the standard MMSN profile. 



cannot at the moment make a strong case against the existence of conditions necessary for GI. 
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Fig. 2. — Protoplanetary core masses at the end of oligarchic 
growth for the standard core accretion model and the 'hybrid', 
migration-enhanced model. M^ng is greater than the predicted 
solid bodymasses for Uranus and Neptune at roughly their present 
orbits in the latter: all of their growth could conceivably occur 
during oligarchy. 
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Fig. 3. — Formation timescales for IOM0 protoplanetary cores 
for 1km and 100m mean accreted planetesimal sizes. The enhance- 
ment of Sp after particle subdisk GI due to Epstein drag-induced 
migration reduces the core formation timescale by nearly lOx for 
dispersion-dominated planetcsimals. lOm-sized planctcsimals drop 
into the shear-dominated regime for a>12.5AU and dominate pro- 
toplanct accretion at these distances such that Tacc ~ lO^yr to the 
edge of the planetesimal disk if 10% of them arc lOm. Thus, for 
s<100m proto-Ncptune can form on timescales Tacc ~ T"d<ss<pi {the 
lifetime of the solar nebula) or less, even at 25 AU. 



